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ABSTRACT 

We have previously introduced the parameter a as an indicator of stability torn = 2 
nonaxisymmetric modes in rotating, self-gravitating, axisymmetric, gaseous (a < 0.34) 
and stellar (a < 0.25) systems. This parameter can be written as a = (ft/2) 1 / 2 , 
where t = T/\W\, T is the total rotational kinetic energy, W is the total gravitational 
potential energy, and / is a function characteristic of the topology /connectedness and 
the geometric shape of a system. In this paper, we extend the stability criterion to 
nonaxisymmetric equilibrium systems by determining empirically the appropriate form 
of the function / for ellipsoids and elliptical disks and cylinders. 

For oblate-like ellipsoidal systems, we find that 

= 2y^f r _ j^sin-ie^/e 2 ) ! 
rf L F(sin _1 e, rj 2 /e 2 )^ 

where e is the meridional eccentricity, r\ is the equatorial eccentricity, and F and E 
are the incomplete elliptic integrals of the first and second kind, respectively, with 
amplitude sin -1 e and parameter rj 2 /e 2 . For prolate-like ellipsoidal systems, we find 
an analogous expression that reduces to / = in the limiting case of infinite cylinders. 

We test the validity of this extension of the stability indicator a by considering its 
predictions for previously published, gaseous and stellar, nonaxisymmetric models. The 
above formulation and critical values account accurately for the stability properties of 
m = 2 modes in gaseous Riemann S-type ellipsoids (including the Jacobi and Dedekind 
ellipsoids) and elliptical Riemann disks as well as in stellar elliptical Freeman disks and 
cylinders: all these systems are dynamically stable except the stellar elliptical Freeman 
disks that exhibit a relatively small region of m = 2 dynamical instability. A partial 
disagreement in the case of stellar Freeman ellipsoids in maximum rotation may be 
due to that the region of instability has not been previously determined with sufficient 
accuracy. 
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Subject headings: galaxies: evolution - galaxies: structure - hydrodynamics - insta- 
bilities - stars: formation 



1 Introduction 

In a recent paper (Christodoulou, Shlosman, & Tohline 1994; hereafter referred to as 
Paper I), we proposed a new criterion for stability of rotating, self-gravitating, axisymmetric 
systems to m = 2 perturbations. The criterion was formulated in terms of the angular 
momentum rather than the energy content of a system and was expressed by the conditions 
a < 0.25-0.26 and a < 0.34-0.35 for stability of stellar and gaseous systems, respectively. 
The stability indicator a was defined for uniformly rotating systems by 

" = (i-D 

where t = T/\W\ is the ratio of the rotational kinetic energy to the absolute value of the 
gravitational potential energy, x = il/fij is the ratio of the rotation frequency to the Jeans 
frequency introduced by self-gravity, and / is a function that depends on both the geometry 
and the topology of a system. 

The first expression of a in equation (1.1) is inconvenient for applications because the 
ratio x = ^/^j cannot be determined easily in all cases of interest. For example, in dif- 
ferentially rotating systems, it is not known how a weighted value of Q should be obtained. 
In centrally condensed and/or nonaxisymmetric systems, it is not known how a weighted 
value of ilj should be estimated. The second expression of a in equation (1.1) is not plagued 
by these difficulties but the functional form of the term / must be known for applications 
to systems with various geometrical/topological structures. A typical case where / is not 
generally known is that of systems with multiply-connected regions (see e.g. the toroidal 
models and the Toomre-Zang disks in Paper I). 

In Paper I, / was determined for homogeneous, uniformly rotating, oblate spheroidal 
systems as 

f = ^L^, (1.2) 
VI - e 2 sin e 

where e is the meridional eccentricity and 



A 1 {e) = ^— 3 sin^e ^— . (1.3) 

In particular, / = 1 for disks and / = 2/3 for spheres. It was also demonstrated in Paper I 
that the second expression of a in equation (1.1) with / given by equation (1.2) provides an 
accurate stability indicator for various oblate spheroidal and disk-like models with nonuni- 
form density and differential rotation. Furthermore, the parameter a was found to be more 
sensitive than the ratio T/|W| that appears in the stability criterion proposed for stellar 
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systems (T/|W| ^ 0.14) by Ostriker & Peebles (1973; hereafter referred to as OP) and in the 
analogous criterion for gaseous systems (T/|W| < 0.27). 

The success of a as an indicator of stability to m = 2 modes in various axisymmetric 
systems prompted us to attempt its generalization in the case of simply-connected, nonax- 
isymmetric systems. In §2, we describe an empirical extension of equation (1.2) to oblate-like 
ellipsoids. With this new definition of /, a is still given by the second expression in equa- 
tion (1.1) and its critical values for marginal stability remain unchanged. In the first four 
subsections of §3, we discuss the predictions of the parameter a for the stability of sev- 
eral oblate-like gaseous and stellar models (ellipsoids and elliptical disks) that have been 
previously constructed and studied by various researchers. In §3.5, we discuss prolate-like 
geometries and the predictions of a for infinite elliptical cylinders. In §4, we summarize our 
conclusions. 



2 The parameters a and / for Oblate-Like Ellipsoids 
Motivated by the results described in Paper I, we adopt the expression 



as a general definition of the parameter a and we obtain a general expression for the function 
/ valid for oblate-like ellipsoids. This method avoids the use of the equations a = t/x and 
X = (see §1 above). This is convenient since it is not clear what form should be 

adopted for the Jeans frequency flj in nonaxisymmetric systems. 

We notice that the terms Ai(e)/(1 — e 2 ) 1 / 2 and e/ sin -1 e in equation (1.2) have been 
introduced by the gravitational potential $ and the gravitational potential energy W, re- 
spectively (Paper I). For oblate-like ellipsoids with three unequal axes a > b > c rotating 
about the short axis c, Ai is a function of both the meridional eccentricity e = (1 — c 2 /a 2 ) 1//2 
and the equatorial eccentricity r\ = (1 — b 2 /a 2 ) 1 / 2 (Chandrasekhar 1969; hereafter referred 
to as EFE), i.e. 



Me, V ) = 2 ^~^f~^ . [ F(sin -i e> 2 /e 2 } _ E (sin -i e> 2 /e2) ] 



(2.2) 



where F and E are the incomplete elliptic integrals of the first and second kind, respectively, 
and 

Ftsin- 1 e, r/ 2 / e 2 ) ee & Jt= , (2.3) 

• 7o y/l - {rf/e 2 ) sin 2 <j> 

/•sin -1 e ; 

^(sin- 1 e, T] 2 /e 2 ) = j sjl- (r] 2 /e 2 ) sin 2 <j> d<j>. (2.4) 

In definitions (2.3) and (2.4), the upper limit of the integrals sin -1 e is called the amplitude 
while the term ?7 2 /e 2 is called the parameter (Abramowitz & Stegun 1972). 
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In addition, W oc F(sin _1 e, r/ 2 /e 2 )/e in oblate-like ellipsoids (e.g. Freeman 1966b) 
and W oc F(7r/2,r/ 2 ) in elliptical disks (e.g. Freeman 1966c; Weinberg & Tremaine 1983). 
Combining the above results, we assume that equation (1.2) generalizes as 

= Me,v) e 

where Ai(e,T]) is given by equation (2.2). For applications, we adopt the expression 

77^ L F(sin e,r) I /e I ) i 

which is derived by combining equations (2.2) and (2.5). For a given model of known 
t = T/\W\, we determine the stability parameter a from equations (2.1) and (2.6) and, as 
in Paper I, we adopt a < 0.25-0.26 and a < 0.34-0.35 to indicate stability to m = 2 modes 
in stellar and gaseous systems, respectively. 



FIGURE 1. The function f is plotted versus equatorial eccentricity 77 for elliptical disks 
(e = 1) and for Jacobi ellipsoids (EFE). 

As we shall see below, equation (2.6) resolves a well-known discrepancy in the case 
of "needles" where e, r\ — >■ 1. Such objects are stable to m = 2 perturbations (e.g. EFE; 
Tremaine 1976; Weinberg 1983). In the limit e, r\ — >■ 1, rj — >■ 1/2 and thus any criterion based 
on t alone fails to predict stability. On the other hand, equation (2.6) shows that / — >■ 
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in this limit and, thus, equation (2.1) predicts stability to m = 2 modes. This behavior is 
seen in Figure 1 where the function f(rj) is plotted for elliptical disks (e = 1) and for Jacobi 
ellipsoids in which e and r\ are uniquely related to each other (e.g. EFE; Christodoulou, 
Kazanas, Shlosman, & Tohline 1994, hereafter referred to as CKST). 

For axisymmetric systems with rj = 0, the above equations reduce to the expressions 
given in Paper I. For nonaxisymmetric disks, e = 1 and the elliptic integrals above become 
complete with amplitude tt/2 and parameter rj 2 (Abramowitz & Stegun 1972; Weinberg & 
Tremaine 1983). Finally, the corresponding expression for the "weighted" Jeans frequency Vtj 
can be determined from the second equality in equation (1.1) and equation (2.6). However, 
such an expression is of limited use since Qj turns out to be a complicated function of rj 
mainly because of the complexity in the functional form of the kinetic energy T in stellar 
systems. 

Analogous equations valid for prolate-like ellipsoids are discussed in §3.5 below in con- 
junction with infinite elliptical cylinders. 



FIGURE 2. Contours of the ratio t = T/\W\ for gaseous Riemann S-type ellipsoids are 
plotted in the (x, b/a) plane. All these objects are stable to second-harmonic perurbations 
(EFE). In the needle limit b/a ^ 0, t 1/2 and is thus insensitive to the dynamics of stable 
strongly nonaxisymmetric Riemann ellipsoids. 
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3 Previously Studied Nonaxisymmetric Systems 

3.1 Gaseous Riemann S-type Ellipsoids 

These uniformly rotating incompressible models are studied in detail in EFE. A brief 
description has recently been given also by CKST. The ellipsoidal figures are oblate-like in 
shape with axes a > b > c and the rotation takes place about the c axis with frequency Q. 



FIGURE 3. Contours of the parameter a for gaseous Riemann S-type ellipsoids are plotted 
in the (x,b/a) plane. All these objects are stable to second-harmonic perurbations (EFE). 
This result is confirmed by the parameter a since a < 0.34 everywhere. In the needle limit 
b/a — > 0, a — > showing its sensitivity to the dynamics of stable strongly nonaxisymmetric 
Riemann ellipsoids. 

All gaseous Riemann S-type ellipsoids, including the Jacobi and Dedekind ellipsoids, 
are stable to second-harmonic perturbations. This behavior is not captured by the stability 
criterion t = T/\W\ < 0.27 that appears to predict instability in the needle limit e, 77 — >■ 1 
where t — >■ 1/2. The criterion based on the parameter a does not suffer from similar difficul- 
ties. The results are illustrated in Figures 2 and 3 where contour plots of the parameters t 
and a are shown in the (x, b/a) plane. The axis ratio b/a = (1 — r) 2 ) 1 / 2 while the parameter 
x = [ab/(a 2 + b 2 )]((/&), where ( is the vorticity in a frame rotating with frequency Q in 
which the ellipsoidal figures appear to be stationary. 
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There exists only one point of marginal stability in Figures 2 and 3 at b/a = 1, x = 1, 
where t = 0.2738 and a = 0.3410 (see also Paper I and CKST). This point denotes the 
onset of second-harmonic dynamical instability in Maclaurin spheroids. All the remaining 
points represent stable Riemann S-type ellipsoids including the stable Maclaurin spheroids 
(b/a = 1). (In fact, all sequences of Riemann S-type ellipsoids characterized by different 
values of x or (/Q bifurcate from the stable part of the Maclaurin sequence.) The parameter 
a, in the form adopted in §2 above, confirms the stability of all these objects to second- 
harmonic perturbations since a < 0.34 everywhere — even in the region \x\ > 5 not covered 
by Figure 3. 



3.2 Gaseous Riemann Disks 



The structure and secular evolution of these uniformly rotating compressible models 
with equatorial axes a > b have been studied by Weinberg & Tremaine (1983). The com- 
pressible Riemann disks and the incompressible ellipsoids of §3.1 exhibit very similar dy- 
namical properties. In particular, all equilibrium sequences bifurcate from the stable part of 
the corresponding Maclaurin sequence of circular disks and all objects are stable to m = 2 
perturbations (Weinberg 1983). 



FIGURE 4. Contours of the ratio t = T/\W\ for gaseous Riemann disks are plotted in the 
(C/fl,b/a) plane. All these objects are stable to m = 2 perurbations (Weinberg 1983). In 
the needle limit b/a ^ 0, t — » 1/2 and is thus insensitive to the dynamics of stable strongly 
nonaxisymmetric Riemann disks. 
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Despite a coordinate change from x in ellipsoids to (/Q in disks (that is done to avoid 
presenting nearly identical diagrams to Figures 2 and 3 and to depict the region near x = 
Q/Vt = in more detail), the dynamical similarities between Riemann disks and ellipsoids are 
apparent in a comparison of the behavior of t = T/\W\ between Figure 2 above and Figure 
4 that shows a contour plot of t for disks in the ((/Q,b/a) plane. An analogous contour 
plot of t in the (b/a, Q 2 ) plane, separated into two figures, can also be found in Weinberg & 
Tremaine (1983). 



FIGURE 5. Contours of the parameter a for gaseous Riemann disks are plotted in the 
((/Q,b/a) plane. All these objects are stable to m = 2 perurbations (Weinberg 1983). 
This result is conhrmed by the parameter a since a < 0.35 everywhere. In the needle limit 
b/a — > , a — > showing its sensitivity to the dynamics of stable strongly nonaxisymmetric 
Riemann disks. 

In the needle limit rj — >■ 1 (b/a — >■ 0), the ratio t = T/\W\ — >■ 1/2 and the stability 
criterion t < 0.27 would predict again instability. In contrast, the parameter a, shown in 
Figure 5 in the ((/fl, b/a) plane, confirms the stability of all Riemann disks to m = 2 modes. 
There exists only one point of marginal stability in Figures 4 and 5 at b/a = 1, (/Q = 2, 
where t = 1/4 and a = 0.3536. In exact analogy to the ellipsoids of §3.1, this point denotes 
the onset of m = 2 dynamical instability in circular Maclaurin disks (see Paper I for more 
details). All the other points in Figure 5 represent dynamically stable disks since a < 0.35. 
This is true also for the disks lying in the region |C/^| > 5 not covered by Figure 5. 
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3.3 Stellar Elliptical Freeman Disks 

These models were constructed by Freeman (1966c) and their stability to m = 2 modes 
was investigated by Tremaine (1976). Brief discussions of the stability properties of Freeman 
disks can also be found in Hunter (1974) and Fridman & Polyachenko (1984). Following 
Hunter (1974), we plot in Figure 6 contours of the ratio t = T/\W\ in the (x 2 ,&/°0 pl ane - 
Here, a > b, x = ^ is the (mean) rotation frequency of the coordinate frame, and flj 

is the equatorial Jeans frequency that depends on the function Ai(e = 1 , 77) . 

This definition of Q in stellar systems with circulation superimposed to the mean rota- 
tion (Freeman 1966b, c; Hunter 1974) is not equivalent to the definition adopted in §3.1 for 
gaseous systems (see Paper I). The difference is also reflected in the behavior oft in the two 
types of systems (CKST). For example, it explains why the Dedekind sequence bifurcates 
from the Maclaurin sequence at t = in stellar systems but at t=0.125 and i=0.1375 in 
gaseous disks and spheroids, respectively. 



FIGURE 6. Contours of the ratio t = T/\W\ for stellar Freeman disks are plotted in the 
(x 2 ,b/a) plane. A roughly triangular region of m = 2 dynamical instability at the upper 
corner of the diagram is marked approximately by the t = 0.14 contour (Tremaine 1976). 
As in the gaseous nonrotating Maclaurin disk = —2, b/a = 1 in Figure 4), t = for 
the Q-model with b/a = 1 and x 2 = 5/8 (cf. Kalnajs 1972). Unlike in gaseous disks, t — >■ 
aiso in the Dedekind limit x —> implying that different definitions of t and Q are used in 
the two cases (CKST; Paper I). 
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In Figure 6, the segment b/a = 1, 5/8 < x 2 < 1 (0 < T/\W\ < 1/2) describes the 
axisymmetric il-models studied by Kalnajs (1972) and Kalnajs & Athanassoula (1974). 
Note that the dimensionless rotation frequency xk that appears in these works (and in §2.3 
of Paper I) is related to the x used here by xk = x(%X 2 — 5)/3; for details see Tremaine (1976). 
The segment b/a = 1, < x 2 < 5/24 (0 < < 0.1286) describes slowly rotating circular 

stellar disks that are secularly unstable because of slow mass loss (hence also slow angular 
momentum loss). These circular disks evolve toward the "stellar-disk Dedekind sequence" 
X = where they finally relax on the segment 0.4 < b/a < 1 becoming nonaxisymmetric 
and nonrotating (see Hunter 1974 and CKST). In the final state, the persisting azimuthal 
flow is interpreted as circulation (Freeman 1966c), hence it does not contribute to the mean 
rotational kinetic energy resulting in T = and t = 0. 



FIGURE 7. Contours of the parameter a for stellar Freeman disks are plotted in the (x 2 ,b/a) 
plane. The roughly triangular region of m = 2 dynamical instability at the upper corner of 
the diagram (Tremaine 1976) is now marked approximately by the a = 0.25 contour. Notice 
that ol — >■ in the needle limit b/a — >■ 0, in the Dedekind limit x —> 0, and for the Q-model 
with b/a = 1 and x 2 = 5/8 (cf. Kalnajs 1972). 

Tremaine (1976) has discovered that m = 2 dynamical instability appears only in a 
roughly triangular region at the upper right corner in Figure 6. The fi-models of Kalnajs 
(1972) with b/a = 1, 5/6 < x 2 < 1 (0.1286 < T/\W\ < 1/2) mark the upper boundary of the 
unstable region. This region is also bounded approximately by the T/| W| = 0.14 contour and 
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by the segment 0.7296 < b/a < 1, x 2 = 1- The point of marginal stability at b/a = 0.7296, 
X 2 = 1 is in very good agreement with the point a/b = 1.3707 of marginal stability in gaseous 
"disk-like Riemann ellipsoids of type I" (EFE). This coincidence is discussed in detail by 
Tremaine (1976). It occurs because in the maximum-rotation limit X = 1 the properties of 
the two types of systems and the corresponding definitions of Q become identical. 

A contour plot of the parameter a in the (x 2 ,b/a) plane is shown in Figure 7. The 
triangular region of m = 2 dynamical instability at the upper right corner of this plot is now 
bounded approximately by the a = 0.25 contour. The agreement between this value and 
the contour value t = 0.14 in Figure 6 is not just an expected coincidence due to the use of 
equation (2.1) with e=l and f{rj) ~ 1 (see Tables 1 and 2 in Paper I). We have compared 
the variations of the two parameters along the marginal stability line given by Tremaine 
(1976). The variation of a along this line is significantly smaller than that of t. Specifically, 
Tremaine (1976) finds that t varies between 0.1286 (at b/a = 1, x 2 = 5/6) and 0.1446 (at 
b/a = 0.7296, x 2 = !)• These values correspond to a variation in t of about 11%— 12%. We 
find that a varies between 0.2536 and 0.2385 at the corresponding end-points. These values 
produce a variation in a of about 6%. The comparison indicates that, as in the axisymmetric 
systems of Paper I, the parameter a is more sensitive than t as a stability indicator of m = 2 
modes in elliptical disks. 

Furthermore, the contour plot of a (Figure 7) confirms the stability of needles unlike 
the ratio T/|W| of the OP criterion (Figure 6). Notice, in particular, that a — > smoothly 
in the needle limit b/a — >■ as well as in the Dedekind limit x ~~ 0. Finally, there exist 
two points of marginal stability with a = 0.2536 along the line b/a = 1 in Figure 7. The 
corresponding critical values are x 2 = 5/24 and x 2 = 5/6. As was described above and in 
CKST, the segment x 2 < 5/24 with t < 0.1286 denotes the appearance of secular instability 
in circular disks unlike the segment x 2 > 5/6 with t > 0.1286 which denotes the appearance 
of dynamical instability (see also Hunter 1974). 

3.4 Stellar Freeman Ellipsoids 

These models of stellar homogeneous ellipsoids were constructed by Freeman (1966b). 
The Freeman ellipsoids are "balanced" in the sense that the gravitational and centrifugal 
forces are by assumption exactly equal to each other along the major axis. This ensures 
that x = 1, where x is defined as in §3.3 for stellar disks. The stability of these models 
to m = 2 modes is not known with certainty and should be further investigated. Fridman 
& Polyachenko (1984; hereafter referred to as FP) describe an attempt to locate regions of 
dynamical instability in the (b/a,c/a) plane, where a > b,c are the principal axes and c is 
the rotation axis of the ellipsoids. Some of their results are suspect because the marginal 
stability curve is jagged. 

To the extent that balanced Freeman ellipsoids share some common properties with 
X = 1 Freeman disks in the regime c/a << 1, we expect that flattened ellipsoids with 
b/a ~ 1 should be dynamically unstable to m = 2 modes (cf. the unstable disks with 
0.7296 < b/a < 1 in §3.3). This expectation is generally consistent with the results shown 
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in Figure 48 of FP but provides no justification for an additional extension of the unstable 
region found at intermediate values of b/a and c/a. It is this extension that appears to be 
suspect. 



FIGURE 8. Contours of the ratio t = T/\W\ for balanced (x = 1) Freeman ellipsoids 
are plotted in the (b/a, c/a) plane. The boundary of a vertical band of m = 2 dynamical 
instability with b/a > 0.73-0.90 (FP) lies between the t = 0.14 and t = 0.20 contours on the 
right-hand side of the diagram. (Thus, the OP criterion is not very successful in tracking 
this boundary.) An extension of the unstable region to intermediate b/a, c/a values shown 
in FP corresponds to low t-values. If real, this extension cannot be understood in terms of 
the OP stability criterion. For disks (c = 0), the exact marginal stability point is located 
at b/a = 0.7296 where t = 0.1446 (§3.3). Virtually all models with b/a < 0.4 have t > 0.14 
although many of them are stable according to FP and all of them are stable according to 
parameter a (see Figure 9). 

Contour plots of t = T/\W\ and a for Freeman ellipsoids are shown in Figures 8 and 
9, respectively. In the limiting case c = 0, the results reduce to those discussed in §3.3 
for Freeman disks in maximum rotation (x = 1). The critical values of both parameters 
(t ~ 0.14 and a ~ 0.25) track the boundary of a vertical band of m = 2 instability shown in 
FP at b/a > 0.73 — 0.90 (for c/a = — 1, respectively) although a is obviously more accurate. 
None of the parameters predicts an extension of the unstable region to intermediate values 
of b/a and c/a seen in Figure 48 of FP. Finally, t is unable to confirm the stability of models 
with b/a < 0.4 since t > 0.14 over most of this region. 
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FIGURE 9. Contours of the parameter a for balanced (x = lj Freeman ellipsoids are plotted 
in the (b/a, c/a) plane. The boundary of a vertical band ofm = 2 dynamical instability with 
b/a > 0.73 — 0.90 (FP) is marked approximately by the a = 0.25 contour. (The variation of 
ol along the boundary is about 5% as opposed to 28%-38% for t in Figure 8.) All models on 
the left of the a = 0.25 contour are predicted to be stable to m = 2 modes by the parameter 
a but not by the results in FP that show an extension of the unstable region roughly between 
the a = 0.15 and a = 0.18 contours. If real, this extension cannot be understood in terms 
of the ol stability criterion. For disks (c = 0), the exact marginal stability point is located 
at b/a = 0.7296 where a = 0.2385 (§3.3). 

3.5 Stellar Elliptical Freeman Cylinders 

These infinite cylindrical models were constructed by Freeman (1966a) and their stability 
properties were investigated by Nishida & Ishizawa (1977) who found them stable to m = 2 
modes. Hunter (1974) presented a contour plot of a parameter t equivalent to the ratio 
T/IlT 7 !. We have repeated the calculation and have reproduced Hunter's result. The contour 
plot of t is shown in Figure 10 in the (x 2 , c /b) plane, where we now assume that the axes 
a > b > c, a — >■ oo, and a is also the rotation axis. We see that the OP criterion, based on 
Hunter's parameter t, suffers from the usual problem of predicting an m = 2 instability for 
a large subset of models with t > 0.14. 
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FIGURE 10. Contours of Hunter's (1974) parameter t (that is equivalent to T/\W\) for 
stellar Freeman cylinders are plotted in the (x 2 , c/b) plane. For these infinite cylinders we 
assume that the rotation is about the a axis, a — > oo, and the equatorial axes b > c. Although 
t > 0.14 in a substantial part of the diagram, all Freeman cylinders are stable to m = 2 
modes (Nishida k, Ishizawa 1977). 



The formulation of §2 is not valid for prolate-like ellipsoidal and cylindrical models that 
rotate about the longest axis. For prolate-like ellipsoids with axes a > b > c and rotating 
about the a axis (EFE), the relevant function is A 2 (e,r/) instead of Ai(e,rj). This function 
is associated with the longest equatorial axis b and can be written as 

Me, V) = 2 f ( ~ V2) l{ 2 [Eism- 1 e, k 2 ) - ^^F^ e, k 2 ) - -^== ] , (3-1) 
r\ l {e l — r\ l ) L e l 1 — r\ l e^l — ri 21 

where now e = (1 — c 2 /a 2 ) 1//2 and n = (1 — c 2 /6 2 ) 1//2 since a is the longest axis and c is the 
shortest axis. The parameter k 2 in the elliptic integrals of equation (3.1) is defined by 

k " = f h ~ V \ Y ( 3 - 2 ) 
By analogy to the assumptions made in §2, we are led to adopt the expression 
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for the function f(e,rj) in prolate-like ellipsoids. Equation (3.3) is analogous to equation 
(2.5) above. The absence of a term similar to (1 — e 2 ) 1//2 of equation (2.5) is due to the change 
in roles of the axes and is justified by a direct calculation in the case of prolate spheroids 
(see below). 

The limiting case of an infinite elliptical cylinder with a — > oo is of interest here in 
relation to the models shown in Figure 10. In this case, e — >■ 1 implying that the amplitude 
sin -1 e — >■ 7r/2 and the parameter k 2 — > 1. Hence, the elliptic integrals in equation (3.1) 
take the asymptotic values E — > 1 and F — >■ oo but the function A 2 (e = l,r)) remains finite. 
By equation (3.3), this asymptotic behavior leads to / = for all infinite cylinders. Hence, 
a = by equation (2.1), indicating stability to m = 2 modes for all Freeman cylinders 
as well as for all gaseous infinite cylinders. This last conclusion is in agreement with the 
result of Chandrasekhar (1981) that circular incompressible infinite cylinders are stable to 
all nonaxisymmetric perturbations including the m = 2 modes. 

The special case of prolate spheroidal models may be of some interest despite the 
fact that a particular class of objects, the homogeneous incompressible prolate Maclaurin 
spheroids, are not equilibrium figures (Florides & Spyrou 1993). For this type of geometry 
and for uniform rotation, we assume that the axes a > b = c and that a is the rotation axis 
and we follow the procedure outlined in Paper I. We find again that t = \f(VL/Vtj) 2 and 
a = (ft/2) 1 ! 2 but now 

/ = (3-4) 



where, as in EFE, e = (1 — c 2 /a 2 ) 1 ! 2 and 



Me) = ^- 1 —^lnJ 1 -±^. (3.5) 



e 2 e 3 



Equation (3.4) is analogous to equation (1.2) above. It has no meaning for homogeneous 
incompressible prolate spheroids but, like equation (1.2), it may prove useful when prolate 
spheroids — stellar or with nonuniform density and/or differential rotation — are examined 
for stability to m = 2 modes. In the limiting cases of e = (sphere) and e = 1 (infinite 
circular cylinder), we recover the expected asymptotic values / = 2/3 and / = 0, respectively. 



4 Summary 



In §2, we have generalized empirically the criterion derived in Paper I for stability 
of rotating, self-gravitating, axisymmetric, stellar (a < 0.25-0.26) and gaseous (a < 0.34- 
0.35) systems to m = 2 modes. The generalized criterion has the same critical values for 
nonaxisymmetric systems as well. As in the axisymmetric case, the fundamental parameter 
a is determined from the equation 

« = /|/*. 

where t = T/\W\ is the ratio of the rotational kinetic energy to the absolute value of 
the gravitational potential energy and / is a function that describes the topological and 
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geometrical structure of a system. For the oblate-like ellipsoids studied in this paper, / is 
given by the generalized expression 



where e and rj are the meridional and equatorial eccentricities and F and E are the incomplete 
elliptic integrals of the first and second kind with amplitude sin -1 e and parameter if/e 2 
(Abramowitz & Stegun 1972). 

In §3, we have compared the predictions of the above formulation with the known 
stability properties of nonaxisymmetric, gaseous and stellar models that have been studied 
previously by various researchers. The parameter a predicts accurately the stability of 
ellipsoids and elliptical disks and cylinders. Specifically, gaseous Riemann S-type ellipsoids, 
including the well-known Jacobi and Dedekind ellipsoids, and gaseous Riemann disks are 
found to be stable to m = 2 modes (in agreement with the results in EFE and in Weinberg 
1983). Stellar elliptical disks (Freeman 1966c; Hunter 1974) exhibit only a relatively small 
region of dynamical instability which lies at the corner of the circular limit rj = and the 
maximum-rotation limit (in agreement with the results of Tremaine 1976). Finally, stellar 
elliptical cylinders (Freeman 1966a; Hunter 1974) are found to be stable to m = 2 modes (in 
agreement with the results of Nishida & Ishizawa 1977). 

This last prediction of the parameter a extends also to gaseous elliptical cylinders (cf. 
Chandrasekhar 1981 where circular incompressible infinite cylinders are found to be stable to 
nonaxisymmetric modes). The reason is that the prediction derives from the condition / = 
(hence a = 0) that is valid for all types of infinite cylinders. This result was obtained in the 
infinite-cylinder limit from the relevant expressions for the function / in prolate spheroids 
and prolate-like ellipsoids. The formulation has been presented in §3.5. 

Little is known about the stability of "balanced" (i.e. maximally rotating) Freeman 
(1966b) stellar ellipsoids with axes a > b,c (see FP). The stability of these models to m = 2 
perturbations should be investigated further. A band of instability extending to all values 
of the axes ratio c/a is shown in Figure 48 of FP in the region b/a > 0.73 — 0.90. This band 
is generally consistent with the results of Tremaine (1976) and with the prediction of the 
parameter a in the disk limit c — > (both indicating an m = 2 dynamical instability in the 
region b/a > 0.73 at maximum rotation). The instability band is predicted accurately for 
all values of c/a by the parameter a (e.g., at b/a > 0.90 for c/a = 1) but not sufficiently 
well by the OP stability criterion. Also consistent with the disappearance of the instability 
in infinite cylinders, the band width shrinks toward zero as c/a increases. However, an 
additional extension of the unstable region to intermediate values of the axes ratios b/a and 
c/a shown in FP is not understood in the present context (see Figure 9 above) and needs to 
be examined in future work. 



/ 



if 



E(sin l e,rf/e 2 ) 
F(sin _1 e, rj 2 /e 2 ) 



(4.2) 
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